Order induced by dipolar interactions in a geometrically frustrated antiferromagnet 
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We study the classical Heisenberg model for spins on a pyrochlore lattice interacting via long 
range dipole-dipole forces and nearest neighbor exchange. Antiferromagnetic exchange alone is 
known not to induce ordering in this system. We analyze low temperature order resulting from the 
combined interactions, both by using a mean-field approach and by examining the energy cost of 
fluctuations about an ordered state. We discuss behavior as a function of the ratio of the dipolar 
and exchange interaction strengths and find two types of ordered phase. We relate our results to the 
recent experimental work and reproduce and extend the theoretical calculations on the pyrochlore 
compound, Gd 2 Ti 2 7 , by Raju et at, Phys. Rev. B 59, 14489 (1999). 



INTRODUCTION 

Geometrically frustrated magnetic systems have re- 
ceived a great deal of attention indecent years from 
both experimentalists and theorists. El Such systems are 
typically composed of corner- or edge-sharing frustrated 
units, usually triangles or tetrahedra. For example, 
a two-dimensional network of corner-sharing triangles 
forms the kagome lattice while the three-dimensional py- 
rochlore structure is composed of corner-sharing tetra- 
hedra. Magnetic systems with such structure exhibit 
unusual low temperature properties which are not com- 
pletely understood. These materials typically lfejmain dis- 
ordered to a freezing temperature T g <C l^cvH-'-™ Indeed, 
the depression of the transition temperature below that 
expected from the measured value of |#cw| is often used 
to gauge the degree of frustration in a magnetic system. 

The reluctance of geometrically frustrated antiferro- 
magncts to order can be understood by considering clas- 
sical models, many examples of which yield macroscop- 
ically degenerate ground states. For instance, it has 
been shown that the Heisenberg model for spins on a py- 
rochlore lattice interacting via nearest neighbor antifer- 
romagnetic exchange has an extensive number of degrees 
of freedom in the ground stateJjQ The degeneracy of the 
ground state manifold is not lifted in this model by ther- 
mal fluctuations and-the system has no finite temperature 
ordering transition.tffi Several realizations ,o£ pyrochlore 
magnets have been studied experimentally.Hflij 

Classical models which include nearest neighbor ex- 
change interactions alone may not be sufficient for ex- 
plaining all the low temperature properties of some frus- 
trated antifcrromagncts. For instance, the origin of the 
freezing transition must involve additional aspects of the 
system. In real materials, several other types of inter- 
actions are present. Further neighbor, exchange may be 
relevant!] as-,may crystal field effectsE3, magnetdc-dipole 
interactionsa or the effect of quenched disordeiotil. 

In this paper we focus on how the inclusion of dipolar 
forces affects the properties of the pyrochlore antiferro- 
magnet. Because of the ground state degeneracy of this 
system with only nearest neighbor exchange, dipolar in- 



teractions are important in establishing order even if they 
are weak. Moreover, the influence of dipole interactions 
in a pyrochlore antiferromagnet has been probedO in the 
compound Gd2Ti2 0y, which should be well represented 
by an isotropic Heisenberg model. 

The properties of Gd2Ti207 have been studied both 
experimentally and theoretically in a recent paper by 
Raju et at, in which it was shown that Gd2Ti2 07 un- 
dergoes a transition— to long range order at a temper- 
ature of about 1K.E3 Their measurement of the high 
temperature susceptibility gives a negative Curie- Weiss 
constant, &cw — ~ 9.8K.E3 Measurements on a magneti- 
cally dilute sample show a reduction in \9cw\, indicating 
that this vaJuc is predominantly due to antiferromagnetic 
cxchangc.113 The transition temperature is much lower 
than l^cwli indicating that the system is frustrated. The 
theoretical work of these authors involves mean-field cal- 
culations expressed Landau expansion of the free 
energy and taken to quadratic order. They examine the 
ordering instabilities that occur as the temperature is 
lowered. With nearest neighbor exchange and long-range 
dipolar interactions, they find that order parameter fluc- 
tuations on entire branches in g-space along the star of 
the_{lll) direction become unstable simultaneously at 
T C S3 Degeneracies of this kind are often broken by ther- 
mal or quantum fluctuations, a phenomenon known as 
order-by-disordcr.Ej Raju et al. suggest that this mech- 
anismj-inay operate to induce order in the model they 
study.E3 

In the following, we extend this mean-field description 
of the system with only nearest neighbor exchange and 
long range dipolar interactions to find the ordering pat- 
tern below T c . We show that the quartic term in the 
free energy expansion lifts the degeneracy of the critical 
modes. The ordering pattern obtained for the ratio of 
dipolar to exchange interaction strengths appropriate for 
Gd2Ti2 07 is a four-sublattice Neel state. In addition, we 
analyze the low temperature fluctuations away from the 
four sub-lattice ground state. We show that all distor- 
tions have a positive energy cost. Ordering in this model, 
therefore, has an energetic origin and is not an example 
of fluctuation-induced order. 
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I. SPINS ON A SINGLE TETRAHEDRON 

It is instructive first to consider spins on a single tetra- 
hedron. Labeling the spins by the interaction energy 
for the tetrahedron is 

Uint E S.j • Sj 

+ J dd -^ ^[S, • - 3(S, • r tf )(Sj • fy)] (1-1) 

where J ex and Jdd are the relative interaction strengths 
of the exchange and dipolar terms, A — fi 2 /a 3 , a is the 
edge length of the tetrahedron, and /i is the dipolc mo- 
ment. We minimize this energy using a standard numer- 
ical search and find all the possible ground states. These 
are shown in Fig. [l[ In each of the ground state configura- 
tions, the spins are coplanar and are tangent to the sphere 
circumscribing the tetrahedron. It is worthwhile to note 
that these states are also ground states for spins inter- 
acting via nearest neighbor antiferromagnetic exchange 
only. Without dipolar interactions, the condition for a 
configuration to be a ground state is that the vector sum 
of the spins on each tetrahedron is zero, leaving two inter- 
nal degrees of freedom for the configuration.Ertl Dipolar 
interactions fix these two degrees of freedom and those 
arising from the 0(3) symmetry. 

Since the interaction energy for a single tetrahedron 
involves only nearest-neighbor dipolar contributions, one 
might not expect the spin configurations in Fig. |l| to be 
a useful guide to behavior on the full pyrochlore lattice. 
In fact, and somewhat surprisingly, the ground state for 
the full lattice that we find in Sec. Ill turns out to be a 
periodic repetition of that for a single tetrahedron, pro- 
vided the ratio of dipolar to exchange interactions does 
not exceed a critical value. 



II. MEAN-FIELD THEORY FOR INTERACTIONS 
ON THE PYROCHLORE LATTICE 

We consider a system of spins interacting via the 
Hamiltonian 

k = \ E s?J?sf 

with 

t q /3 x x i ^ ix 3(r?j)(rfj) 

J if = O a pd n ,i + — (dap 2 ) ( 2 -l) 

where C = u 2 Jdd/ 'Jex- The exchange interaction is be- 
tween nearest neighbors only, as indicated by S nn . The 
labels I and j refer to sites in the lattice and a and (3 label 
components of the spin or spatial vectors. A mean-field 




FIG. 1. Three degenerate ground states for spins on a sin- 
gle tetrahedron. The spins in each configuration are parallel 
to certain edges of the tetrahedron. These edges and the spins 
parallel to them are drawn with the same type of line, either 
bold or dashed. There are three other ground states, obtained 
from these by reversing all spins. 



treatment for Hamiltonians of this kind has been devel- 
oped by Reimers, Berlinsky and Shi.Q For completeness, 
we briefly summarize their approach. Order parameters 
for the system are defined by 

Tr( ft S, ; ) = B, (2.2) 

where pi is the local density matrix and is constrained by 
Tr(pi) = 1. Expanding the free energy in a power series 
in E$i yields 

F = const. 

+\ E Bf{iT5 a0 8 l3 +Jtf)B^ 

+ ^Tj2(B?B?)(B?BP) + (D(B e ). (2.3) 

The next step is to diagonalize the quadratic term. In 
systems with n atoms per unit cell, it is convenient to 
divide the site label i into two parts: I labels the unit 
cell and a labels spins within the unit cell. Making the 
Fourier transforms 

5f Q =E B q Qelq ' ri ( 2 - 4 ) 

q 

— J_ jaab/3 -iq-r im In r\ 

J lm.ab - N 2^ J q 6 ^- a > 

q 

and substituting these expressions into the free energy, 
one arrives at the expression for the quadratic part of 
/ = F/N, the free energy per unit cell, used in Ref. [l2|: 

/ (2) = ^EC {iT5 ab 5 a() + J^) B%. (2.6) 

ab,a/3 q 

Consider diagonalization of the 'in x 3n matrix, J^ ab ^ . 
Denote the eigenvalues by Aq and eigenvectors by U q for 
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2 = 1... 3n. Expanding the order parameters in the basis 
of eigenvectors 



one obtains 



(2) _ 



EE 

q % 



3T + A' 



(2.7) 



(2.8) 



V 4-> 



To find the minimum of the quadratic term, it is neces- 
sary to determine the minimum A q for all i and q. 

Specializing to the pyrochlore lattice, we follow previ- 
ous workers in using a rhombohedral unit cell contain- 
ing four magnetic ions (3n = 12)Jj'E3 Basis vectors are 
! ( j + k) , | (i + k) , and | (i + j) , where a is the edge length 
of the cubic unit cell. Also, the four magnetic ions within 
the unit cell are located at (x,y,z) = (0,0,0), (0, j, 
(i,0,i), and (i,i,0). 

The instability of the paramagnetic phase is analyzed 
by considering the sign of (3T + A q ). In the paramagnetic 
phase, all are positive. At the transition, the smallest 
becomes negative. Using the ratio Jdd/Jex — 0.2, Raju 
et al. showed that such an approach does not completely 
determine the ordering pattern.113 At T c , a star of modes 
simultaneously becomes unstable. For the same value 
of Jdd/Jex , we recover these results. By contrast, for 
values of Jdd/Jex > J c — 5.7, we find that a discrete 
set of isolated modes (related by the lattice symmetry) 
become unstable. 

To illustrate this point, we plot in Figure |^ the mini- 
mum eigenvalue of the J q matrix as a function of q for 
various values-of J = Jdd/Jex ranging from that used 
by Raju et alt3 (based on the measured values of \9cw\ 
and the lattice spacing and the calculated value of the 
magnetic moment of Gd 3+ ) to a much larger value. For 
J < J c , the minimum eigenvalue is independent of q 
along the (111) direction. For J > J c , there are isolated 
minima located close to (but not at) q = 0. These are 
the individual modes that become unstable at T c in the 
case J > J c . 

We now investigate the ordering pattern below T c for 
J < J c . The degeneracy of the soft modes is lifted 
when we include contributions to the free energy from the 
fourth order term. At temperatures below T c , the order 
parameters acquire finite magnitude and one is faced with 
the problem of simultaneously minimizing the quadratic 
and quartic parts of the free energy. Consider first the 
quartic term in isolation. In terms of the real space order 
parameters, Bf a , this is 



9T 
20 



Ed« jDOt d/5 r>/5 
la La la la 

la, a/3 



(2.9) 



Its value depends on the magnitude and direction of 
B; a . For fixed magnitude of T,i aa (B" a ) 2 , the quartic 
term is minimized by a state with all |Bj a | equal. For- 
tunately, and apparently fortuitously, a state satisfying 
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FIG. 2. Minimal eigenvalues along the (111) line in q-space 
for various values of the dipole-dipole interaction strength. 
Curve a corresponds to Jdd/Jex = 0.2 as in Ref. 12. Curves b 
through / have J dd /J ex = 2.92, 5.63, 8.35, 11.06, and 13.78, 
respectively.tJ 




FIG. 3. An illustration of the eigenvectors for the soft 
modes along the star of the (111) directions in q-space. The 
vector q is perpendicular to the base of the tetrahedron and 
is shown coming out of the page. The three spins are copla- 
nar, have equal magnitude, and are tangent to the cirle which 
circumscribes the base of the tetrahedron. The spin on the 
fourth site has zero magnitude. 



this condition can be constructed which also minimizes 
the quadratic term. In detail, we proceed as follows. 
The eigenvector associated with one of the modes that 
becomes unstable at T c is illustrated in Fig.||. Taking 
the most general combination of this eigenvector and its 
three partners, and imposing the condition that the val- 
ues of |B; a | are the same at four sites of a tetrahedron, 
we generate the configurations of Fig.^ and no others. 
Tiling the lattice with these configurations, we obtain 
only states with ordering wavevector q = 0. We con- 
clude that the ordering pattern for T just below T c is as 
shown in Fig. [l| 
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FIG. 4. A projection of our q = ordering pattern onto 
the xy-p\ane. The spins are coplanar and form a 4-sublattice 
Neel state. 



III. STATIC DISTORTIONS AWAY FROM THE 
ORDERED STATE 

Finally, we investigate whether a ground state with this 
ordering pattern is stable against thermal fluctuations 
at low temperatures. To describe this calculation using 
the notation introduced above, we impose the constraint 
|Bj | = 1 for every I and a. We denote the components of 
B; a in the ground state with the ordering pattern shown 
in Figure [| by Bi aa . We expand 



Blaa — B laa 



6B laa -±B? aa \SB la \ 



0(<5B 4 ). (3.1) 



Since <5B; a (when small) is orthogonal to B/ Q , it has only 
two independent components. We expand the energy to 
quadratic order in the distortion variables obtaining 



mb{3 



lmab 7 a{3 

l 

2 



~(\6B la \ 2 + \5B mb \ 2 )B? aa Bi b p) (3.2) 



where Eq is the energy of the ground state. We then 
diagonalize the quadratic term, which involves an 8 x 8 
matrix for each wavevector, q. We find the minimum 
eigenvalues of this matrix as a function of q. As ex- 
pected, we recover a flat branch of zero-energy modes 
when J dd — 0.0. □ For J dd > 0.0, we find that all fluctu- 
ations away from an ordered state are associated with a 
positive energy cost (even in the long wavelength limit, 
since dipolar interactions break global rotational symme- 
try). As an illustration, we plot the minimum at each q 
along the (001) direction in q-space for various values of 
the dipole-dipole interaction in Figure |. The energetic 
barriers around our proposed ground state at T = mean 
that there is not a degenerate, connected manifold. The 
ordering obtained is due to energetic selection and does 
not occur via an order-by-disorder mechanism. 
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FIG. 5. The minimum eigenvalues for distortions away 
from the ground state in the (001) direction in q-space. The 
various curves correspond to different values of the relative 
strength of the dipole-dipole interaction, Jdd/Jex- Curves A 
through E have J dd /J ex = 0.0, 0.04, 0.08, 0.12, and 0.16, 
respectively. Curve F has Jdd/Jex = 0.2 as in Ref. 12. 



SUMMARY 

We have considered the influence of dipolar interac- 
tions in a nearest neighbor Heiscnberg model of a geo- 
metrically frustrated systeMi reproducing and extending 
earlier work by Raju et. alxl In particular, we determine 
the ordering pattern, using a mean field treatment and 
by examining the stability of an ordered state to fluctu- 
ations. We believe that the ordering we find should be 
that associated with the transition which_ia observed ex- 
perimentally at around IK in Gd2Ti 2 07.li3 This ordered 
state, shown in Figure ||, is a four sub-lattice Neel state. 

It will be interesting to compare these theoretical 
predictions with the results of neutcpn scattering from 
160 Gd2Ti2O7, currently in progress .113 
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